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We discuss the appareance of ordered structures in assemblies of magnetic particles. The phe- 
nomenon occurs when dipolar interactions and the thermal motion of the particles compete, and is 
mediated by screening and excluded volume effects. It is observed irrespective of the dimensional- 
ity of the system and the resulting structures, which may be regular or fractal, indicate that new 
ordered phases may emerge in these systems when dipolar interactions play a significant role. 
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I. INTRODUCTION 

The collective behavior of assemblies of dipolar mag- 
netic particles is governed by the influence of several com- 
peting mechanisms and presents a very rich and com- 
plex phenomenology whose ultimate interest is the 
identification of the emerging structures and the sub- 
sequent analysis of its implications in the macroscopic 
properties of the system. When the particles are dis- 
persed in an aqueous, magnetically inactive, phase, the 
formation of structures results — as occurs in other soft- 
condensed matter systems — from the interplay of dipolar 
interactions, excluded volume effects, and thermal mo- 
tion. Unlike assemblies of non-magnetic particles, they 
may exhibit ordered phases emerging when dipolar inter- 
actions are dominant, and disordered phases when those 
interactions are inhibited by Brownian motion. In the 
cases of electro- and magneto-rheological dispersions, the 
particles are heavy enough to impede Brownian effects, 
and therefore chaining is the most likely event. Contrar- 
ily, for finer magnetic particles, as monodomain particles 
or ferrofiuids, thermal motion becomes more important, 
opening the possibility to the analysis of the interplay 
with dipolar interactions, and the subsequent formation 
of structures. Disconnected results along this lines have 
been previously reported for different systems [^j3| . 

Our purpose in this paper is to discuss an appar- 
ent general phenomenon occurring in these systems: the 
emergence of order when magnetic interactions and ther- 
mal motion compete. Ordered and disordered states 
would then be accessible upon variation of temperature, 
and consequently would correspond to different struc- 
tural phases of the system. 

The paper is distributed as follows. In Section II we 
study the problem of aggregation of dipolar particles in 
d = 2 and d = 3, whereas in Sec. Ill we focus on a prob- 
lem of different nature, namely the adsorption of particles 
onto a substrate of dimension d = 1 and d — 2. Finally, 
in the last Section we discuss the common features of 
the transition to more ordered structures when dipolar 
interactions play a significant role. 



II. CLUSTER AGGREGATION WITH DIPOLAR 
INTERACTIONS 

The first system we consider is the cluster aggregation 
of particles that experience dipolar interactions. Cluster 
aggregation, and in general fractal growth phenomena 
are very active fields, specially in which regards 
the physics of colloids. In the context of computer mod- 
els, the most notable are the diffusion-limited aggrega- 
tion (DLA) model Q and the cluster-cluster aggregation 
model 0,11 , which effectively describe the fractal struc- 
ture of real colloidal aggregates . In this Section, we 
describe the extension of the standard DLA model taking 
into account fully anisotropic dipolar interparticle inter- 
actions. 



A. Description of the algorithm 

The general framework of the models we are consider- 
ing is the sequential addition of dipolar particles of diam- 
eter a, with a moment /x rigidly attached, to an assembly 
S of rigid dipoles located at fixed positions. The incom- 
ing particles perform a random walk under the influence 
of the long-range interactions exerted by the dipoles in 
S. The effect of these interactions is encoded in the fol- 
lowing procedure: Consider that, at a certain time step, 
the structure S is formed by N rigid dipoles /Xj, located 
at positions Ri. The energy of an incoming particle of 
magnetic moment fi at position r, is then given by pO[ 
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where = r — Ri is the relative distance between par- 
ticles. In the next time step, the dipole is moved to a 
new position r* = r + Sr, where Sr is a small increment 
oriented along a randomly chosen direction. The energy 
in the new point r* is £* , and the total variation in the 
energy due to the movement is AS = £*—£. If Af < 0, 
the movement is accepted with probability 1; otherwise, 
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it is accepted with probability exp(— A£a^//i^T,.), where 
Tr = a^kBT/ fj? is a dimensionless reduced temperature, 
relating the intensity of the dipolar interactions /it to the 
actual temperature T of the system. In case of not being 
accepted, the trial position r* is discarded and a new one 
generated. 

With the previous algorithm, we have taken into ac- 
count the translations dynamics of the particles. For the 
rotational dynamics of the dipoles, we apply the follow- 
ing prescription: The initial orientation of the dipoles ijP 
is assigned at random. After every accepted movement, 
the moment of the random walker is oriented along the 
direction of the total magnetic force on its position. Af- 
ter the dipole becomes finally attached to the assembly, 
its moment undergoes a final relaxation, and its direction 
does not change any more. Indeed, this fact assumes that 
the relaxation time for the orientation of a particle is very 
short compared with the characteristic time scale of the 
movement of its center of mass. 

The simulation of the particle-cluster aggregation pro- 
cess starts with a particle (the seed) located at the origin 
of coordinates, and bearing a randomly oriented moment 
IXq. The following particles are released from a random 
position on a sphere of radius TZi^ centered on the seed. 
The particles have an initial moment fj.^ assigned at ran- 
dom. Each particle undergoes a Brownian motion, ac- 
cording to the given prescription, until it cither contacts 
the cluster or moves away a distance larger than T^out- 
In this case, the particle is removed and a new one is 
launched. In the present simulations, we have used the 
values 7?.out = 27?. and Ttin = T^out — 5a, where TZ is 
the radius of the cluster, measured in units of the par- 
ticle diameter a. The particles stick to the cluster when 
they overlap one or more particles previously aggregated. 
The new dipole is then attached in the point of its last 
trajectory where it is first tangent to the surface of the 
cluster. 



B. Numerical results in d = 3 

Simulations of cluster aggregation of dipolar particles 
in d = 3 are extremely time-consuming. Therefore, we 
report only results for the limits = and Tr — oo. 
The case Tr — Q corresponds to infinitely strong dipolar 
interactions; the steps during the random walk are only 
accepted if they lower the total dipolar energy. For Tr = 
oo all movements are accepted and we must recover the 
standard DLA model 

In Figure Q we show two typical clusters grown at 
Tr ~ and = cxd in a box of size L — 150a. The growth 
was stopped when the radius of the clusters reached a 
predetermined distance to the boundary of the box. In 
the case Tr = 0, the cluster was composed by iV ~ 1500, 
whereas for = oo we have TV ~ 10000. A first conclu- 
sion can be immediately drawn: At high temperatures, 
the clusters are more compact, having an roughly spheri- 




a) 



b) 



FIG. 1. Projection onto the x — y plane of two typical dipo- 
lar clusters in d = 3, grown on a cubic box of size L = 150a. 
(a) Tr = oo, D/ = 2.47 ± 0.03. (b) Tr = 0, _D/ = 1.37 ± 0.03 
(see text). 



cal shape. When decreasing the temperature, the average 
density of the clusters diminishes; at = 0, the clusters 
have an essentially digitated structure of very low den- 
sity. 

A more quantitative estimate of this concept can be 
obtained by measuring the fractal dimension -D/ of the 
clusters The value of D/ is determined from a log- 
log plot of the radius of gyration Rg{N) as a function 
of the number of particles in the cluster N, using the 
relationship B 
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The value computed for the aggregation at = oo is 
Df = 2.47 ± 0.03, in good agreement, within error bars, 
with the best estimate for DLA in rf = 3, namely D^^^ = 
2.495 ±0.005 0. On the other hand, for the simulations 
at zero temperature, we obtain a much smaller value, 
Df = 1.37 ±0.03. 
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These strong variations in the structure of the clus- 
ter with temperature can be understood through a very 
simple and intuitive argument. The aggregation process 
depends on the competition between dipolar forces and 
thermal agitation. Dipolar interactions favor the forma- 
tion of chains of aligned dipoles, thus minimizing the 
dipolar energy, while the thermal motion tends to ran- 
domize the process and produce a higher degree of ramifi- 
cation. This branching process is mediated by the screen- 
ing of the inner regions of the cluster by the most external 
branches. The chances for a random walker to surpass 
the external active region and reach the core of the cluster 
are very small, and this effect induces a high branching 
ratio and a large fractal dimension. The balance between 
these two effects, branching and growth, can be measured 
by estimating the relative probabilities of growth, pg, of 
an already existing branch and its splitting, ps, creat- 
ing thus a new branch. Assuming that the branches are 
composed by relaxed (parallel) dipoles, pointing along 
the axis of the branch, the addition of a new particle on 
the tip of an existing branch will increase the energy of 
the aggregate by an amount A£g = — (we only 

consider the interaction of the new particle with its near- 
est neighbor). We can associate to this growth event a 
relative probability given by a Boltzmann factor, and set 
Pg = exp(— AfgC^/^^Tr) — exp(2/rr). On the other 
hand, in order to minimize the energy, an incoming par- 
ticle, inducing a branching event, will set its dipole anti- 
parallel to its nearest-neighbor in the existing branch, 
thus increasing energy by ASg — —[i^jc?. The corre- 
sponding splitting probability is then ps — exp(l/rr). 
We have therefore 



^ = exp(l/TO. 

Vs 
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At high temperature, pg ^ Ps', the growth and splitting 
of a given branch are equally likely events. The effect 
of the dipolar interactions is overcome by the thermal 
disorder and we must recover the original DLA with no 
interactions. At low temperatures, the splitting probabil- 
ity is vanishingly small compared with pg , and one should 
expect to observe in this limit very fingered aggregates, 
with little branching, and a fractal dimension closer to 1, 
as is indeed the case in Fig. |l|. 

Considerations about the fractal dimension allows also 
to draw conclusions on the macroscopic aggregation of 
individual clusters. At low concentrations of dipolar par- 
ticles, it is reasonable to assume that particles first form 
clusters, which on their hand coalesce later to form larger 
structures. If the mutual interactions between clusters 
are negligible (due to screening of the interactions in a 
large reservoir of clusters with random total moment) 
the coalescence of clusters is mediated by Brownian mo- 
tion. At low temperatures, the cluster dimension is low. 
Let and D^p be the fractal dimension of any two 
clusters in the reservoir (of course, bmce 
d'^P + d'^P < 3, the clusters are transparent, have little 
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FIG. 2. Fractal dimension Df for dipolar clusters in d = 2 
as a function of the reduced temperature Tr- 



mutual interactions and behave as essentially indepen- 
dent objects At high temperatures, however, the 
fractal dimension is larger and fulfills the condition of 
mutual opacity + > 3 Assuming that Df 
is an increasing function of Tr (as is the case in d = 2, see 
next section), we predict the existence of a critical tem- 



perature Tr for which 
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Df ' = 3, setting the thresh- 



old above which cluster coalescence is possible, leading to 
the formation more complex structures. 



C. Numerical results in d = 2 

Simulations in c? = 2 are much less costly, and afford 
the possibility of truly exploring the effect of tempera- 
ture on the cluster structure. Moreover, simulations on 
a plane are also realistic and interesting on their own, 
because of their relevance in understanding phase transi- 
tions in Langmuir monolayers p3[ | . Far from equilibrium, 
the condensed phase grows at the expense of the liquid 
phase, forming clusters of different shapes. The phospho- 
lipids that make up the Langmuir monolayer experience 
repulsive dipolar interactions, which must play a major 
role in determining the morphology of the condensed ag- 
gregates. 

In Figure |^ we plot the fractal dimension D f for dipo- 
lar clusters in = 2 as a function of the reduced tem- 
perature Tr- This figure confirms the arguments given in 
Sec. [IB. At high temperatures the fractal dimension is 
large {Df — 1.74 ±0.03 at Tr = 10) and compatible with 
the best estimates for free DLA {D^^^ = 1.715 ± 0.004, 
pl|). At low temperatures, on the other hand, the frac- 
tal dimension is much smaller, and closer to 1 {Df — 
1.13 ± 0.01 at Tr = 0). That is, the lower the tempera- 
ture, the stronger the screening effects of the branches, 
leading to an essentially digitated structure composed 
by fingers with very few side branches. At this respect, 
a significant difference with the case 3 is that now 



3 



0.40 



0.36 



^ 0.32 



0.28 



0.24 





o T, 


= 




AT, 


= 10-2 






= 10 








1 


1 


1 1 



JC/4 



71/2 



37t/4 



FIG. 3. Angular correlation function g{8) in dipolar clus- 
ters in ci = 2 for different values of the reduced temperature 



clusters are always mutually opaque, D^j'' +Dy' > 2, for 
all values of Tr. 

Another remarkable feature of this system is the or- 
dering induced by the dipolar interactions. A convenient 
way to measure this order is the angular correlation func- 
tion g{9) of the dipoles composing the cluster, defined as 
the fraction of all pairs of dipoles, irrespective of the po- 
sition, forming a relative angle 9. This correlation func- 
tion is shown in Fig. ^ for different temperatures. At 
high temperatures we observe that the system is com- 
pletely disordered, with all orientations between dipoles 
being equally probable. More interestingly, at low tem- 
peratures we observe the emergence of an ordered phase, 
in which dipoles are parallel with high probability, while 
suppressing the anti-parallel (more energy costly) config- 
urations. That is, there is an order-disorder transition 
(reflected in the angular correlation function g) between 
an ordered phase at the low-temperature limit (in the ab- 
sence of any thermal disorder) in which the orientation of 
the dipoles is strongly correlated, and a disordered phase 
in the high-temperature limit, in which dipoles show an 
uncorrelated behavior, imposed by the intrinsic fractal 
geometry of the clusters. 

The structures that we have shown in this section have 
been observed experimentally [^4|jl^. In the latter ex- 
periments, ferromagnetic particles of different sizes are 
transferred to a Langmuir monolayer at different pres- 
sures. The main observed feature is that the larger the 
diameter of the particles, the more digitated the emerging 
cluster. These experimental results can be interpreted by 
using our argument which is able to predict the shape of 
the clusters in terms of the reduced temperature. That 
equation indicates that, under those circumstances, the 
probability of a cluster to grow is much higher than the 
probability for splitting. Conversely, when the diameter 
of the particles is small, both probabilities equilibrate, 
resulting in denser circular compact aggregates, also in 



accordance with experimental results. 



III. CHAINING OF DIPOLAR PARTICLES ONTO 
A SUBSTRATE 

In this Section we present results for a model of adsorp- 
tion in which the adsorbing particles experience dipolar 
interactions. The irreversible adsorption of colloidal par- 
ticles onto a solid surface is a very relevant subject, with 
plenty of practical and technological applications ||l^,0 ■ 
The basic models proposed so far in an attempt to un- 
derstand the physics of adsorption are the random se- 
quential adsorption model ||l^ and the ballistic model 
(BM) jl^. The BM model is a good approximation for 
the adsorption of large colloidal particles. In this model, 
particles arrive at the surface following straight vertical 
trajectories. When an incoming particle fails to reach 
the substrate directly, it is allowed to roll down over the 
previously adsorbed ones, until it reaches an equilibrium 
position. Particles eventually resting on the surface are 
adsorbed. Otherwise, they are rejected. 

These models, and their subsequent generalizations, 
consider only very short range — usually hard-core — in- 
teractions. Here we will illustrate the profound effect of 
interparticle dipolar interactions in the adsorbed phase, 
for the case of ballistic adsorption. 



A. Description of the algorithm 

We consider the irreversible adsorption of dipolar par- 
ticles that diffuse in a semi-infinite d + 1-dimensional 
space and adsorb upon contact onto a boundary d- 
dimensional hyperplane (the substrate). In this case, the 
assembly of particles S under consideration is the phase 
formed by all the particles adsorbed at given time, lo- 
cated at positions Ri and carrying a rigid moment fi^. 
Since we are considering an irreversible process, the ad- 
sorbed particles are fixed and cannot diffuse on top of the 
substrate. The algorithm used to simulated the random 
walk of the incoming particles is exactly as described in 
section II A, differing of course in the boundary condi- 
tions. 

Simulations start with a clean substrate. Particles are 
released sequentially at random positions over a certain 
initial height h-m above the substrate, and with an initial 
random orientation /Xq. The particles perform a random 
walk until they either hit the substrate and become ad- 
sorbed or move away a distance larger than hont- In 
this case, the particle is removed and a new one re- 
leased. We have chosen the parameters /lin = 15a and 
hout — 20a, which are reasonable given the power-law 
decay of dipolar interactions (jl|) . The rotational dynam- 
ics of the dipoles is the same as in the case of cluster 
aggregation, namely, dipoles relaxing at every accepted 
step. When a particle reaches a position very close to the 
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FIG. 4. Saturated configurations in d = 2, system size 
L = 10a. (a) Free BM, Tr = oo. (b) With dipolar inter- 
actions, at Tr = 0. The particles extend out of the box to 
show the effect of periodic boundary conditions. 

surface (less than one particle diameter a) , it is attached 
to the substrate following the BM rules. This procedure 
assumes therefore the presence of a strong short-ranged 
interaction between substrate and particles, which over- 
comes dipolar interactions at very short distances. The 
newly adsorbed particles experience one last relaxation, 
and their moment does not change any more afterwards. 
The adsorbing surface has linear dimension L and we im- 
pose periodic boundary conditions in all directions but in 
the perpendicular to the plane. The system sizes consid- 
ered in this work range between L ~ 150a and L ~ 25Ga. 

B. Numerical results in d = 2 

Again, numerical simulations in d = 2 are very time- 
consuming and permit only to consider very small system 
sizes. In Figure ^ we have represented two typical sat- 
urated configurations, corresponding to the adsorption 
according to the rules of standard BM {Tr — oo) and 
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FIG. 5. Chain-length density function P{£) for different 
values of the reduced temperature Tr, for dipolar adsorption 
in d = 1. System size L = 250a. 

in the presence of dipolar interactions (T^ = 0). Inspec- 
tion of these figures allows to conclude that the packing 
of the substrate is higher in presence of interactions. A 
numerical measure of the packing is given by the jam- 
ming limit 9oo, defined as the maximum fraction of the 
surface covered by particles at the saturation limit. Our 
simulations provide a value Ooo = 0.622 ± 0.003 for the 
dipolar adsorption at = 0, whereas the standard BM 
{Tr = oo) corresponds to 9^ = 0.61048 ± 0.00013 [|o). 
This increase in the packing of the saturated surface can 
be traced back to the attractive interactions between the 
adsorbed phase and the incoming particles. In the case of 
standard free BM, the structure of the adsorbed phase is 
due exclusively to excluded volume effects: the already 
attached particles impede the adsorption of new parti- 
cles at random positions, and allow only the occupation 
of empty spaces. This is essentially a random, disor- 
dering, effect. The presence of dipolar interactions, at- 
tracting the random walker to the surface, diminishes 
the excluded volume effects and leads to configurations 
of higher packing. 

C. Numerical results in d = 1 

By restricting the adsorption onto a line, we can ex- 
plore more thoroughly the interplay between interactions 
and the disordering excluded volume effects. As we can 
conclude from examination of Fig. ^, the reason of the 
higher coverage is the larger tendency of dipolar parti- 
cles to form connected structures at low temperatures. 
In d = 1 these connected structures are identified with 
chains. We therefore propose as a measure of order the 
chain-length density function P{(.), defined as the aver- 
age number of chains of length la, per unit length of 
substrate. From this function we can compute the re- 
maining interesting properties of the adsorbed phase, like 
the jamming limit 
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FIG. 6. Jamming limit 6*00 as a function of the reduced 
temperature Tr, for dipolar adsorption in d = 1. 
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In Figure ^ we have represented the chain-length den- 
sity function for dipolar adsorption at = and 
Tr = 00 (free BM). At high temperature we observe a 
very fast (super-exponential) decay of P{t), indicative 
of a lack of any characteristic length; disorder has sup- 
pressed the formation of correlated chains in this case. At 
low temperatures, however, P{(.) decays exponentially, 
P{t} ~ exp(— with a well-defined correlation length 
(or average chain size) ^ = 4.73 ± 0.02. Chains have now 
the possibility of growing very large, and a single chain 
is seen to cover more than 15% of the whole substrate. 

Finally, we have plotted in Fig. |^ the jamming limit 
as a function of the reduced temperature T^. At high 
temperature [Tr = 10) we recover the theoretical value 
6»BM ^ 0.8079 ± 0.0003 |l9). On the other hand, at 
low temperatures the higher packing due to the chain- 
ing of the particles leads to a higher coverage, namely 
^'oo = 0.855 ± 0.002, which differs from the free case in 
about 6%, a notable variation in the context of adsorp- 
tion phenomena. Note, moreover, the strong resemblance 
of Figs. ^ and |^, which depict analogous magnitudes in 
both models presented. 



IV. CONCLUSIONS 

The results presented in this paper render the following 
conclusions. The overall behavior of assemblies of mag- 
netic particles is the result of the competition of factors 
of different nature: dipolar interactions, thermal disor- 
der (Brownian motion) and screening (for aggregation) 
or excluded volume effects (for adsorption). 

The concept of reduced temperature enables us to as- 
certain the global form of the clusters. A more precise 
analysis of the geometry of the structures, able to account 
for their local structure, needs however the use of fractal 
concepts, and in particular of fractal dimension. In the 



cases we have discussed, magnetic particles may assemble 
basically into structures of effective Df — 1, correspond- 
ing to the chaining process on a substrate mediated by 
excluded volume interactions, or into structures with val- 
ues of Df comprised in the interval ]1,2[, ranging from 
digitated to spherical clusters. 

In the case of dipolar cluster aggregation, we have ob- 
served the existence of transitions between an ordered, 
or quasi- ordered, phase emerging at low values of the 
reduced temperature (and therefore comprising the low 
temperature and large particle size cases) and a disor- 
dered phase for high values of the reduced temperature. 
The low temperature phase is characterized by the exis- 
tence of long-range correlations between dipoles. In the 
d = 2 case, our predicted results have been corroborated 
in experiments of aggregation performed on Langmuir 
monolayers. This technique has revealed to be a very 
useful tool in the characterization of the resulting two- 
dimensional structures. 

Finally, it is worth pointing out that the existence of 
order at the mesoscopic level may have implications in 
the behavior of macroscopic quantities, such as the mag- 
netization of the system. 
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